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An  Integrated  Experimental  and  Computational  Study  of  Relocalization  in 
7r-Conjugated  Combustion  Intermediates 
Andrew  L.  Cooksy 

San  Diego  State  University  Department  of  Chemistry  and  Biochemistry 

1  Statement  of  Problem  Studied 

Relocalization  isomers  are  configurational  isomers  connected  by  a  vibronic  pathway  that  redistributes 
electron  density,  altering  orbital  hybridizations  and  molecular  geometry  without  changing  the  sequence 
of  chemical  bonds  (Fig.  1).  Appearing  in  hydrocarbon  radicals  and  carbenes,  relocalization  confers  a 
structural  ambiguity  to  these  reaction  intermediates  that  introduces  a  largely  unstudied  complexity  to 
the  chemical  pathways  of  combustion  [1,  2]. 


Figure  1:  Canonical  forms  of  molecules  investigated  in  this  work. 


The  goal  of  this  work  has  been  to  probe  the  dynamics  of  relocalization  coordinates  in  hydrocarbon 
and  oxyhydrocarbon  free  radicals  by  ab  initio  electronic  structure  and  vibrational  calculations  and 
by  hot  band/combination  band  infrared  vibrational  spectroscopy.  We  sought  to  compute  vibrational 
transition  energies  and  intensities,  vibrationally  averaged  geometries,  and  other  dynamical  properties 
for  HC3O  [3]  and  other  radicals  relevant  to  hydrocarbon  combustion  [2,  4],  employing  coupled  cluster 
calculations  and  correlation-consistent  basis  sets  to  generate  accurate,  multidimensional  potential  energy 
surfaces.  From  these  surfaces,  the  coupled  vibrational  Schrodinger  equation  is  solved  using  a  package 
developed  by  our  group  and  based  on  the  finite  element  method.  Furthermore,  experimental  spectra 
can  be  measured  by  adapting  a  tunable  infrared  spectrometer  to  operate  at  liquid  nitrogen  temperature. 

Because  these  intermediates  are  formed  early  in  hydrocarbon  combustion,  their  ability  to  behave  like 
two  chemically  distinct  structures  may  have  a  dramatic  influence  on  current  kinetic  models.  These  stud¬ 
ies  therefore  directly  address  ARO  research  areas  7.1  (Experimental  Physical  Chemistry,  experiments 
and  calculations  that  enable  modeling  of  the  time  dependent  processes  of  ignition  and  combustion)  and 
7.6  (Theoretical  Chemistry:  predictive  computational  methods  for  chemical  processes,  e.g.  combustion). 

2  Results  Summary 

2.1  Computational  Quantum  Mechanics 

Using  new  computational  resources  provided  by  this  grant,  we  have  constructed  the  first  potential 
energy  surfaces  for  relocalization  that  parametrize  the  system  along  all  of  the  local  coordinates  involved 
in  the  motion.  The  first  surface  we  evaluated  is  for  the  HC3O  radical,  using  coupled  cluster  methods  and 
triple-zeta  quality  basis  sets,  evaluating  the  energy  at  about  250  points  as  a  function  of  the  HCC  and 
CCO  bending  angles,  and  along  the  CCC  asymmetric  stretching  coordinate.  Figure  2  shows  the  high 
degree  of  anharmonicity  in  the  potential,  including  the  strong  angular-radial  coupling.  These  features 
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Figure  2:  The  HC3O  potential  energy  surface  as  a  function  of  HCC  and  CCO  bending  angles  at  selected 
values  of  R  =  R( C1C2)  —  -R(CiC2). 
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arise  from  the  low-energy  isomerization  pathway  open  to  these  radicals,  which  we  believe  may  lead  to 
early  branching  of  the  combustion  reaction  pathways  of  unsaturated  hydrocarbon  fuels.  An  analogous 
series  of  calculations  has  also  been  carried  out  for  the  OglU  cyclooctatetraenyl  radical,  formed  by  the 
reaction  between  acetylene  and  benzyl  radical  (Fig.  3).  These  surfaces  will  be  improved  by  the  use  of 
new  correlation-consistent  basis  sets  prior  to  publication  in  peer-reviewed  journals.  The  manuscripts 
for  this  work  will  be  submitted  to  the  Journal  of  Chemical  Physics  and  Journal  of  Physical  Chemistry 
A. 

Using  these  surfaces,  we  then  employed  our  own  software,  FENMvib  [5],  to  extract  the  wavefunc- 
tions  and  energies  of  the  vibrational  states  corresponding  to  excitation  of  the  isomerization  coordinate. 
Support  of  this  project  has  led  us  to  extend  the  capabilities  of  FEMvib  to  automatically  interpolate  and 
extrapolate  the  grid-based  potential  energy  surface  by  means  of  distributed  approximating  functionals 
[6].  This  is  a  crucial  step  towards  providing  a  web-based  platform  to  allow  other  researchers  to  solve 
the  vibrational  Schrodinger  equation  in  a  black-box  fashion,  providing  only  their  potential  energy  grid 
and  corresponding  geometries.  Efforts  in  this  direction  continue  in  our  lab. 


Figure  3:  The  CsfU  potential  energy  surface  as  a  function  of  CCC  angle  and  framework  torsion  angle 
at  selected  values  of  the  bond  length  R  which  changes  formally  between  single  and  double  bond  during 
relocalization. 


2.2  Rovibrational  Spectroscopy 

Our  solutions  to  the  vibrational  Schrodinger  equation  reveal  a  dramatic  effect  of  the  flat  surface  on  the 
bending  frequencies  of  these  reactive  intermediates.  In  HC3O,  the  putative  CCH  bending  excitation 
frequency  appears  at  680  cm-1  and  the  CCO  bend  at  345  cm-1.  Furthermore,  extremely  strong  coupling 
of  the  bending  and  stretching  coordinates  leads  to  modes  such  as  the  one  at  680  cm-1  in  which  the 
wavefunction  node  is  along  the  HCC  bending  coordinate  at  high  R  (as  defined  in  Fig.  2)  but  along 
the  CCO  coordinate  at  low  R.  These  low-frequency  bending  modes  and  strong  coupling  enhances  the 
likelihood  of  our  observing  hot  band  and  combination  band  spectra  in  our  spectrometer. 

Support  of  our  related  experimental  work  has  upgraded  an  existing  spectrometer  (Fig.  4)  to  allow 
us  to  scan  for  new  rotationally  resolved  infrared  spectra  of  these  radicals.  We  have  been  scanning  the 
2000  cm^1  region  for  carbon  stretching  transitions  of  HC3O,  using  acetylene  and  carbon  monoxide  to 
make  the  radical  in  an  electric  discharge,  according  to  methods  determined  in  previous  work.  Although 
the  high  rotational  partition  function  has  precluded  observation  of  these  asymmetric  radicals  in  the 
past,  the  liquid  nitrogen  cooling  funded  by  this  grant  should  bring  the  transition  intensities  within  the 
detectable  range  of  our  experiment. 

Upgrades  of  the  spectrometer  are  nearly  complete,  and  we  will  continue  scanning  first  for  HC3O 
transitions,  now  with  improved  estimates  of  the  combination  band  transition  frequencies.  Subsequent 
searches  will  be  carried  out  for  transitions  of  the  acroleinyl  radical  C3H3O,  based  on  forthcoming  analysis 
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Figure  4:  Schematic  of  the  tunable  diode  laser  spectrometer. 


of  the  vibrational  eigenfunctions  and  eigenvalues  of  that  system. 

2.3  Roaming  Free  Radical  Mechanism  in  HC3O 

In  addition  to  the  work  described  above,  the  grant  supported  our  extension  of  the  work  on  HC3O  to 
investigate  this  molecule’s  possible  relevance  to  the  roaming  radical  problem  so  elegantly  investigated 
by  Suits  and  coworkers  [7].  The  HC3O  radical  is  a  chemical  analogue  of  the  formyl  radical,  with  an 
additional  ethynyl  linking  the  oxygen  to  the  terminal  hydrogen.  However,  the  added  chain  length  allows 
HC3O  to  exist  in  either  a  vinylic  form  (bent  CCH,  similar  to  formyl)  or  an  acetylenic  form  (straight 
CCH),  with  the  acetylenic  form  more  stable.  We  wished  to  investigate  the  possibility  that  the  roaming 
hydrogen  extraction  channel  energetics  were  significantly  affected  by  the  relocalization  coordinate.  Our 
work  finds  that  this  is  not  the  case.  Instead,  incipient  removal  of  the  attached  H  atom  by  the  roaming  H 
atom  is  most  stable  along  a  vinylic  geometry,  in  spite  of  the  substantially  different  equilibrium  geometry. 
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